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Abstract 

Let \P] be the points of a Poisson process on and F a probability 
distribution with support on the non-negative integers. Models are for¬ 
mulated for generating translation invariant random graphs with vertex 
set \P\ and iid vertex degrees with distribution F, and the length of the 
edges is analyzed. The main result is that finite mean for the total edge 
length per vertex is possible if and only if F has finite moment of order 
{d + \)/d. 
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1 Introduction 

Consider a Poisson process V on with intensity 1 and write [P] for the point 
set of the process. Furthermore, take a probability distribution F with support 
on N. How should one do to obtain a translation invariant random graph with 
vertex set [V] and degree distribution F7 Which properties does the resulting 
configuration have? These are the questions that will be considered, and partly 
answered, in this paper. 

The problem of generating random graphs with prescribed degree distribu¬ 
tion has been extensively studied the last few years. It is primarily motivated 
by the increasing use of random graphs as models for complex networks; see 
e.g. Durrett (2006). Such networks typically exhibit non-Poissonian distri¬ 
butions for the degrees, predominantly various types of power-laws, and it is 
hence interesting to develop graph models that can produce this type of degree 
distributions (indeed, the standard Erdos-Renyi graph gives Poisson degrees 
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in the limit). One example of a graph model for obtaining a prescribed de¬ 
gree distribution is the so called configuration model - see Molloy and Reed 
(1995,1998) - where, given n vertices, each vertex is independently assigned 
a random number of stubs according to the prescribed distribution and the 
stubs are then paired randomly to create edges. See also Bollobas et al. (2001), 
Chung and Lu (2001:1,2) and Bollobas et al. (2006) for other work in the area. 

Most existing models for obtaining graphs with given degree distribution 
do not take spatial aspects into account, that is, there is no metric defined 
on the vertex set. In Deijfen and Meester (2006) however, a model is intro¬ 
duced where the vertex set is taken to be Z. The vertices are equipped with 
stubs according to a prescribed distribution and these stubs are then paired 
according to a certain rule. This model, which turns out to be closely related 
to random walks, is shown to lead to well-defined configurations, but the edge 
length has infinite mean. In Deijfen and Jonasson (2006), a different model 
(on Z) is formulated and shown to give finite mean for the edges. This work 
is generalized in Jonasson (2007) to more general (deterministic) vertex struc¬ 
tures, in particular it is shown there that there exists a translation invariant 
graph on Z'^ with finite total edge length per vertex if and only if the degree 
distribution has finite mean of order {d + l)/ii. 

Here we will take the vertex set to be the random point set \P] of the 
Poisson process V. Each vertex is independently assigned a number of stubs 
according to the desired degree distribution F and we then ask for a way of 
connecting the stubs to form edges. The restrictions are that the procedure 
should be translation invariant and the resulting graph is not allowed to con¬ 
tain self-loops and multiple edges, that is, an edge cannot run from a vertex 
to itself and each pair of vertices can be connected by at most one edge. A 
few different ways of doing this will be discussed, and we will show that the 
condition for the possibility of finite mean for the total edge length of a vertex 
is the same as on Z'^. 

As described above, the model can (at least for d = 1, 2) be seen as a way 
of introducing geography in a complex network model. However, it can also 
be viewed as a generalization of the problem treated in Holroyd et al. (2008), 
where the matching distance in various types of matchings of the points of a 
Poisson process is analyzed. Indeed, this would correspond to taking F = 1 
in our setup, that is, to equip each vertex with one single stub. We will make 
use of results from Holroyd et al. at several occasions. 

To formally describe the quantities we will work with, let (P, ^) be a 
marked Poisson process, where the ground process P is a homogeneous Pois¬ 
son process on with rate 1 and the marks are iid with distribution F; see 
Daley and Vere-Jones (2002, Section 6.4). We will often refer to the points 
in [P] as vertices, and we think of the marks as representing an assignment 
of stubs to the vertices. Let A4 be a translation invariant scheme for pairing 
these stubs so that no two stubs at the same vertex are paired and at most 
one pairing is made between stubs at two given vertices. Furthermore, let 
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be the Palm version of {V,i,M.) with respect to V and write 
P* and E* for the associated probability law and expectation operator respec¬ 
tively, that is, P* describes the conditional law of given that there 

is a point at the origin, with the mark process and the pairing scheme taken 
as stationary background; see Kallenberg (1997, Chapter 11) for details. Note 
that, since the Palm version of a homogeneous Poisson process has the same 
distribution as the process itself with an added point at the origin, we have 
that [P*] = [P] U {0}. Let T denote the total length in M* of all edges at the 
origin vertex. We will show the following result for T. 

Theorem 1.1 There exists a translation invariant pairing seheme withW\T] < 
oo if and only if F has finite moment of order {d+ l)/d. 

Let D denote the degree of the origin. The only-if direction of the theorem 
is fairly easy and follows from investigating the radius of the smallest ball 
around the origin of that contains D points in V*. Indeed, since multiple 
pairing to the same vertex is not allowed, T is bounded from below by the sum 
of the distances to the vertices in this ball, which, as it turns out, cannot have 
finite mean if is infinite. As for the if-direction, we will describe 

a concrete pairing scheme with the property that E*[T] < oo. Here, results 
from Jonasson (2007) and Holroyd et al. (2008) will be useful. 

Before proceeding, we remark that the distribution of T can be described 
in an alternative way: For r € [0, oo], let H{r) denote the expected number of 
points in \P] D [0,1]'^ whose total edge length does not exceed r. Since V has 
intensity 1 it is not hard to see that 77 is a distribution function, and, in fact, 
the distribution of T is given by H; see Holroyd et al. (2008, Section 2) for a 
more careful description of this connection. 

The rest of the paper is organized so that the only-if direction of Theorem 
11.11 is proved in Section 2, a possible pairing schemes is described in Section 
3 and a refinement of this scheme that proves the if-part of Theorem 11.11 is 
formulated in Section 4. Section 5 contains an outline of possible further work. 

2 Only-if 

Recall that D denotes the degree of the origin vertex in V*. We will show 
that E*[r] < oo is impossible if = oo. To this end, let B{r) be the 

d-dimensional ball with radius r centered at the origin and write V*{r) for the 
number of P*-points in B{r), excluding the origin, that is, 

V*{r) = caidiV* nB{r)] - 1, 

where card[-] denotes set cardinality. Furthermore, for n > 1, define Rn to be 
the radius of the smallest ball centered at the origin that contains at least n 
points of P* not counting the point at the origin, that is, 

Rn = inf{r : V*{r) > n}. 
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The idea is that, with a probability that can be bounded from below, at 
least half of the D stubs at the origin will be connected to stubs at vertices 
whose distance to the origin is at least On this event we have that 

T > DRd/{2-3-^/‘^) and, as we will see, the mean of Rn behaves like for 
large n. Combining this will give that T must have infinite mean if 
has so. 

Consider a Poisson process V\ on with intensity A and - in analogy with 
the above definitions - let V{r) = caidlVxriB{r)] and Rn = inf{r : V{r) > n}. 
The mean of Rn is characterized in the following lemma. 


Lemma 2.1 For n > 1, we have that 


lE[^n] 


_c_ A r(t + i) 

xi/d z^r(/t + i)’ 

k=0 


( 1 ) 


where C = C{d) is a constant. In particular, £[!?„] > C for all n > I 

and some constant C = C'{d). 


Proof. Since the volume of B{r) is cr^, where c = c{d) is a constant, we 
have that V (r) is Poisson distributed with parameter Acr*^. Together with the 
definition of Rn, this yields that 


^{Rn >r)= P(P(r) <n) = Y^ 

k=0 


(Acr-^)^ Acr'i 
kl 


Hence 

roo ^ 1 poo 

> r)dr = V TT / (Acr'^)*^e^"^ dr, 

Jo Jo 


and, by variable substitution. 



r(fc + ^) 

(i(Ac)^/'^ 


Since k\ = r(A: + 1), this establishes ([T|) with C = c~^^^d~^. As for the second 
claim of the lemma, just note that, since r(A: + ^)/r(A: + 1) is decreasing in k, 
we have that 

v r(fc + i) ^ r(n + i) r(n + i) 

^jik + i)- ■r(n + i) r(n) ’ 

and, by Sterlings formula, r(n + ^)/r(n) ~ as n —oo. □ 

With Lemma [Q at hand it is not hard to prove that finiteness of E*[T] 
requires finite moment of order (d + l)/d for D. 
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Proof of onlv-if direction of Theorem fl. 11 Write E*[T] = E* [E*[r|Z), 
and, conditional on D and Rd, define 


A = 



< 


D - 1 
2 


that is, A is the event that at least half of the D — \ non-origin P*-points in 
the interior of B{R£)) are located at distance at least from the origin. 

Since V* has the same distribution as a Poisson process with an added point 
at the origin, it follows from properties of the Poisson process that the D — 1 
non-origin points are uniformly distributed in B{Ro), implying that 


E* 


V* {^ 
I 3Vrf 


D, Rd 


= {D-1)- 


yoI[B{3-^/^Rd)] _D-1 
voI[B{Rd)] ~ 3 ^ 


where vol[-] denotes volume. Combining this with Markov’s inequality, we can 
bound 


and hence 


E*[T\D,Rd] > 


{D-l)/2 
E*[T\A,D,Rd] 


On the event A, at least {D — l)/2 of the D stubs at the origin must be 
connected to vertices whose distance to the origin is at least implying 

that 


Hence, using Lemma l2.11 we get that 


> 


D -1 Rd 
- 2 3Vrf 

get that 


6 • 3V<^ 

6 • 3V<^ 

6 • 3^/^ 

1)Rd] 

1)E*[Rd\D]] 


which proves that < oo is indeed necessary for E*[T] < oo. □ 


3 Repeated stable matching 

In this section we make a first attempt to formulate a pairing scheme. The 
algorithm is based on so called stable matchings - see e.g. Gale and Shapely 
(1962) - obtained by iteratively matching mutually closest points, and turns 
out to give finite mean for T in d > 2 if P has bounded support. For d = 1 
an alternative scheme is described. 
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For a general translation invariant homogeneous point process TZ with finite 
intensity, one algorithm for matching its points is as follows; First consider 
all pairs {x, y} C [7^] of mutually closest points - that is, all pairs {x, y} 
such that X is the closest point of y and y is the closest point of x - and 
match them to each other. Then remove these pairs, and apply the same 
procedure to the remaining points of [7^]. Repeat recursively. If TZ is such 
that almost surely [TZ] is non-equidistant and has no descending chains, then 
this algorithm can be shown to yield an almost surely perfect matching, that 
is, a matching where each point is matched to exactly one other point. Here, 
the point set [TZ] is said to be non-equidistant if there are no distinct points 
x,y,u,v G [7^] with |x — y| = [u — u|, and a descending chains is an infinite 
sequence {x*} C \TZ\ such that \xi — Xi_i| is strictly decreasing. Furthermore, 
the obtained perfect matching is the unique stable matching in the sense of 
Gale and Shapely (1962). 

Now consider our marked Poisson process, where each point x G \P] has 
a random number Dx ^ F of stubs attached to it. For each point x, number 
the stubs 1,..., Dx and say that stub number i belongs to level i in the stub 
configuration. Furthermore, for i > 1, define \P]i = {x G \P] : Dx > i}- 
To connect the stubs, first match the points in \P]i to each other using the 
stable matching algorithm described above, and pair the stubs on the first level 
among themselves according to this matching. Then take a uniformly chosen 
2-coloring of the points in \P ]2 such that points that are matched to each 
other in the matching of \P]i get different colors. This gives rise to two sets of 
points [ 7 ^) 2,1 and [7^]2,2 with different colors. Match the points in [ 7 ^) 2,1 to each 
other using the stable matching algorithm and do likewise with the points in 
\P]‘ 2 , 2 - Connect the stubs on the second level accordingly. Then continue in 
this way. In general, the points in \P]i are colored with a uniformly chosen 
i-coloring such that points that have been matched in the previous steps are 
assigned different colors, giving rise to point sets [7^]j,i,..., [P]i,i whose points 
are then matched among themselves using the stable matching algorithm. 
Since the distributions of the point sets {['P]i,j} fulfill the conditions for the 
stable matching algorithm to yield a perfect matching, this procedure indeed 
provides a well-defined pairing of the stubs, and the coloring prevents multiple 
edges between vertices. We refer to this scheme as repeated stable matching 
with coloring (RSMC). 

As for the total edge length per vertex with RSMC we have the following 
result. 

Proposition 3.1 For RSMC with degree distribution F, we have that 

(a) W[T] = 00 for any degree distribution F in d = 1. 

(b) E*[T] < 00 for F with bounded support in d >2. 

Proof. Part (a) is a direct consequence of Theorem 5.2(i) in Holroyd et al. 
(2008), where the expected matching distance for the stable matching of a one- 
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dimensional Poisson process is shown to be infinite. As for (b), write D = Dq 
and let Li denote the length of the edge created by the hth stub at the origin 
(if D < i, then Lj := 0). Also write pi = ¥{D = i) and pj_|_ = ¥{D > i). By a 
straightforward generalization of Theorem 5.2(ii) in Holroyd et al. (2008), we 
have that 


P* (Li > r\D >i)< C—r-^ (2) 

Pi+ 

for some constant C = C{d). Indeed, the point sets are out¬ 

comes of translation invariant processes with intensity Pi+ji, satisfying the 
conditions for the existence of a well-defined stable matching, and, as pointed 
out by Holroyd et ah, the proof of Theorem 5.2(ii) applies to any such process. 
It follows that E*[Lj|D > i] < Ci/piJ^ for some constant C = C{d) and hence, 
with u = max{i : pi > 0}, we have 


U 

E*[T]=E[Y,Li 

i=l 


u u 

^^Pi+E*[Li\D > i] < C i < oo. 

i=l i=l 


□ 


Remark. Note that the proof applies also if the marking of the Poisson points 
with random degrees is not independent, but only translation invariant, that 
is, the probability law of the marked process is invariant under shifts of M'^. 
Indeed, Theorem 5.2(ii) can be applied to conclude ([2]) also in such a situation. 
□ 


Clearly there is no reason to believe that RSMC is optimal to create short 
edges. Firstly, for a Poisson process, the stable matching is not optimal in 
terms of the matching distance - see Holroyd et al. (2008) - indicating that 
RSMC is not optimal within each level. Secondly, RSMC suffers from the 
obvious drawback that stubs on different levels are not connected to each 
other. Indeed, a stub on level i must go at least to the nearest other vertex 
with degree larger than or equal to i to be connected. By Lemma 12.11 the 
expected distance to such a vertex is implying that 


OO / n \ 

^*[T]>^(^Cp-^nPn. 

n=l Vi=l / 

This lower bound is infinite for instance for a power law distribution with 
sufficiently small exponent t > 1 , that is, for a distribution with pn = Cn~'^. 
Indeed, then pn+ = so that pi+ < for i = n/2,... ,n, 

and we get 


E*[T]>cY,^ 

n=l 




oo 

n=l 
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which is infinite for r < {2d — 1) /{d — 1). In d = 2 hence E*[r] = cxd for r < 3 
with RSMC, while Theorem 11.11 stipulates the existence of a pairing scheme 
with E*[r] < oo for r > 5/2. 

Before describing such a pairing scheme, we mention that, in d = 1, RSMC 
can be modified to give finite mean for degree distributions with bounded 
support as follows: As before, let [V]i = {x € [P] : > i}- To connect 

the stubs on level i, pick one of the two possible matchings of the points of 
[Vji satisfying that, for any two matched points x, y with x < y, the interval 
{x,y) does not contain any points of [V]i. Call a point that is matched with 
a point to the right (left) right-oriented (left-oriented) and connect the stubs 
so that the stub at a right-oriented point x is paired with the stub at the hth 
left-oriented point located to the right of x. This is clearly stationary, and, 
since [V]i C [V]i-i, no multiple edges are created. We refer to this procedure 
as shifted adjacent matching (SAM). 

Proposition 3.2 Assume that F has bounded support and let u = max{pj : 
Pi > 0}. For SAM, we have that E*[T] = In particular, E*[T] < oo if and 
only if F has bounded support. 

Proof. Recall that Li is the length of the edge created by the f:th stub 
at the origin {Li := 0 if D < i). With SAM, the stub at level i at the 
origin (if such a stub exists) will be connected to the {2i — l):th point in 
\P]i to the right or left of the origin with probability 1/2 respectively. The 
distance from the origin to the {2i — l):th point to the right with degree at 
least z is a sum of a NegBin(2z — l,Pi+) number of Exp(l)-variables. Hence 
E*[Lj|Z? > i] = {2i — l)/pi+ and the proposition follows from a calculation 
analogous to the one in the proof of Proposition 13.11 □ 

Remark. SAM for a stub configuration on is analyzed in Deijfen and 
Jonasson (2006: Section 2), where it is shown that T has finite mean also for 
stationary degrees {Di} with bounded support. □ 

4 Finite mean is possible when < oo 

We now describe a scheme that gives finite mean for the total edge length 
per vertex as soon as F has finite mean of order {d + f)/d. Roughly, the 
model is a truncated version of RSMC, where vertices with very high degree 
are connected nearby (instead of having to go to other vertices with equally 
high degree) and the remaining stubs are then connected according to RSMC. 
The model is designed to make it possible to exploit results from Z'^ and the 
procedure for connecting the high-degree vertices will involve associating the 
Poisson points with their nearest vertex in a uniform translation of Z'^. 

The following stepwise algorithm for connecting vertices with large degree 
to vertices with small degree in an iid stub conhguration {D^} on Z'^ was 


formulated in Jonasson (2007, Section 3.2). It is a generalization of a similar 
algorithm for Z from Deijfen and Jonasson (2006) inspired by the ’’stable mar¬ 
riage of Poisson and Lebesgue“ from Hoffman et al. (2006), and it is applicable 
also to more general graphs. For a large integer m, let D'^ = Dzl{Dz > m} 
and say that z is high if D' > 0, otherwise z is called low. First, the posi¬ 
tions of the vertices of Z'^ are disturbed by moving each vertex independently a 
Unif(0,0.1)-distributed distance along a randomly chosen incident edge (this is 
just to make the vertex set non-equidistant). Then, in the first step, every high 
vertex z claims its nearest low neighbors, and a low vertex that is claimed 
by at least one high vertex is connected to the nearest high vertex that has 
claimed it. Let ^^^(l) be the number of remaining stubs of the high vertex z 
after this has been done. In the second step, each high vertex z with Dz{l) > 0 
claims its Dz{l) nearest low vertices that have not yet been connected to any 
high vertex, and each claimed low vertex is connected to the nearest high 
vertex that has claimed it. This is then repeated recursively. In Jonasson 
(2007) it is shown that, if m is chosen large enough, this procedure leads to a 
well-defined configuration with E[T] < oo as soon as < oo, where 

T denotes the total edge length (in the Z'^-metric) of the origin. 

Now return to the Poisson setting. For x G M'^, write Ux for the unit cube 
centered at x. Pick a point xq uniformly in the origin cube Uq and let Z'^(xo) 
be a translation of Z'^ obtained by moving the origin to xq. The stubs at a 
Poisson point x G [P] are said to be associated with the point z G Z‘^(xo) 
such that X G Uz- Write Dz for the number of stubs associated with a point 
z G Z'’*(xo) and Nz for the number of Poisson points in Uz- Then Nz ~ 
Po(I) and Dz = Di where {Di} are the iid marks of the Poisson 

points in Uz- The variables {Dz} are also iid and can be thought of as a stub 
configuration on Z'’*(xo). We number the stubs associated with each vertex 
z G Z‘^(xo) randomly in some way, for instance by first numbering the Poisson 
points in Uz randomly 1 ,..., Nz, and then consecutively picking one stub from 
each Poisson point according to that ordering until all stubs are numbered. 
Furthermore, in order to keep track of which Poisson point a certain stub 
originates from, the stubs are labelled with the position of their Poisson point. 
The stubs are now connected as follows: 

1. Pick m large and, as in Jonasson (2007), say that a vertex z G Z‘^(xo) 
is high (low) if Dz > m {Dz < m). Match the stubs of the high vertices 
of Z‘^(xo) with stubs at the low vertices of Z'^(xo) according to the al¬ 
gorithm from Jonasson described above, using the stubs at each vertex 
in numerical order. When two stubs are matched, we create an edge 
between the Poisson points that they originate from. Since there will be 
no multiple matchings between the same two vertices in Z'^(xo), this will 
not give multiple edges between the Poisson points either. 

2. After step I, in d = 1, the unconnected stubs - that is, the stubs at the 
low vertices of Z(xo) that have not been matched with stubs from the 
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high vertices - are connected with SAM. In d = 2, we drop the associa¬ 
tion of the stubs with the vertices of Z‘^(xo) and take the unconnected 
stubs back to their Poisson points. These stubs are then connected with 
RSMC. 

This pairing scheme is clearly translation invariant and will give finite mean 
for T if < oo: li Di has finite moment of order {d + \)/d, then 

has so as well. Hence, if the point z G Z‘^(xo) to which the stubs at the origin 
vertex of \P*] are associated is high, then it follows from Jonasson (2007) that 
T has finite mean. If the point is low, then T has finite mean by Proposition 
IQ and 13.21 and the following remarks. Indeed, the remaining stubs after step 
1 induce a translation invariant assignment of degrees to the Poisson points - 
or, in d = 1, to the points of Z(xo) - where the degrees are bounded by m. 

5 Further work 

We have formulated a necessary and sufficient condition for the existence of 
a random graph on a spatial Poisson process with finite expected total edge 
length per vertex. There are several related problems that remain to investi¬ 
gate. 

Generalization of stable matching. A natural model for pairing the stubs 
is the following: Consider the marked Poisson process and, at time 0, start 
growing a number of balls from each Poisson point linearly in time, the number 
of balls that start growing from a given point being given by the mark of the 
point. When the collection of balls from two points meet, one ball from each 
point is annihilated, and an edge between the two points is created. The 
remaining balls continue growing (this is to avoid multiple edges). When 
D = 1, so that each stub has exactly one stub attached to it, this gives rise to 
the unique stable matching of the points; see Holroyd and Peres (2003, Section 
4). It remains to analyze the algorithm for other degree distributions, which 
seems to be more complicated. 

Percolation. Apart from the edge length, other features of the configura¬ 
tions arising from different pairing schemes could also be investigated. One 
such feature is the component size, that is, the number of vertices in a given 
component of the graph. Will the resulting edge configuration percolate in 
the sense that it contains an unbounded component? This question is not 
relevant when D < 1 - indeed, the configuration will then consist only of 
isolated edges, implying that the component size is at most 2 - but arises 
naturally for other degree distributions. How does the answer depend on the 
degree distribution? For a given degree distribution, does the answer depend 
on the pairing scheme? Is it always possible to achieve percolation by taking 
d sufficiently large? 
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Independent Poisson processes. Instead of a single Poisson process, the 
vertex set could be generated by two independent Poisson processes, repre¬ 
senting two different types of vertices. This is related to matchings of points 
of two independent Poisson processes considered in Holroyd et al. (2008). We 
look for ways of obtaining a graph with edges running only between differ¬ 
ent types of vertices and with prescribed degree distributions for both vertex 
types. Can this be done if the degree distributions are different for the two 
vertex types? If so, how different are the degree distributions allowed to be? 
Which properties do the resulting configurations have? 


Acknowledgement I thank Ronald Meester for discussions on the problem 
during a stay at VU Amsterdam spring 2008. 
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